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Abstract. RR Lyrae stars in the Milky Way are good tracers to study the kinematic behaviour and spatial 
distribution of older stellar populations. A recently established well documented sample of 217 RR Lyr stars with 
V < 12.5 mag, for which accurate distances and radial velocities as well as proper motions from the Hipparcos 
and Tycho-2 catalogues are available, has been used to reinvestigate these structural parameters. The kinematic 
parameters allowed to calculate the orbits of the stars. Nearly 1/3 of the stars of our sample has orbits staying near 
the Milky Way plane. Of the 217 stars, 163 have halo-like orbits fulfilling one of the following criteria: Q < 100 
kms"'^, orbit eccentricity > 0.4, and normalized maximum orbital ^-distance > 0.45. Of these stars roughly half 
have retrograde orbits. The z-distance probability distribution of this sample shows scale heights of 1.3 ± 0.1 kpc 
for the disk component and 4.6 ± 0.3 kpc for the halo component. With our orbit statistics method we found 
a (vertical) spatial distribution which, out to 2; = 20 kpc, is similar to that found with other methods. This 
distribution is also compatible with the ones found for blue (HBA and sdB) halo stars. The circular velocity Q, 
the orbit eccentricity, orbit z-extent and [Fe/H] are employed to look for possible correlations. If any, it is that 
the metal poor stars with [Fe/H] < 1.0 have a wide symmetric distribution about = 0, thus for this subsample 
on average a motion independent of disk rotation. We conclude that the Milky Way possesses a halo component 
of old and metal poor stars with a scale height of 4-5 kpc having random orbits. The presence in our sample of 
a few metal poor stars (thus part of the halo population) with thin disk-like orbits is statistically not surprising. 
The midplane density ratio of halo to disk stars is found to be 0.16, a value very dependent on proper sample 
statistics. 
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K*" 1. Introduction 



Studying the spatial distribution of stellar populations in 
the Milky Way requires well defined tracer objects. Among 
the large variety of stars, those of horizontal-branch (HB) 
nature are well suited. HB stars are core-He burners after 
the red-giant phase whose envelopes have little (< 0.1 Mq) 
to perhaps 0.7 Mq, leading to a location on the HB from 
the very blue to red (sdB, HBB, HBA, RR Lyr, RHB 
stars). They had an initial mass on the main sequence 
of < 2.5 which means that some of the RRLyrs are 
only about 1 Gyr old and therefore belong to the younger 
disk population. In particular the blue HB stars and RR 
Lyrae variables are easy to find. For the RR Lyrae stars 
(RR Lyrs) the mean absolute magnitude is well calibrated 
using data from the Hipparcos (de Boer et al. 1997a; 
Fernley et al. 1998). The distances of blue HB stars are 
easily determined utilising their well defined My - For the 
blue HB stars, models based on the evolutionary nature 
allow to derive Tcs and log g from their spectra, thus lead- 
ing to a fairly accurate My as well. Red HB stars are less 
easy to identify, but the Hipparcos data allow to iso- 



late nearby ones from the parallax based CMD (Kaempf 
et al. 1995). 

As a result of many years of work on the nature 
and distribution of sdB stars (Heber 1986; Moehler 
et al. 1990; Theissen et al. 1993; Saffer et al. 1994, 
Villeneuve et al. 1995) the idea emerged to include radial 
velocity and proper motion data allowing to calculate the 
kinematical properties of these stars, including the calcu- 
lation of orbits (following Odenkirchen & Brosche 1992). 
Several such studies have been carried out thus far. 

Using sdB stars, de Boer et al. (1997b) derived from 
the orbit statistics of 41 stars a z-distribution with a scale 
height of ~ 1 kpc plus a hint of the existence of a more ex- 
tended population. The existence of this latter population 
was substantiated by Altmann et al. (2004) in a sample 
of 114 stars, showing a halo group with a z-distribution 
which one can characterise by a scale height of ~ 7 kpc. 
Recently, using the CMD identification of red HB stars, 
Kaempf et al. (2005) presented a similar study of the RHB 
star distribution. They find a disk-like distribution with a 
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z-scale height of perhaps up to 0.8 kpc and an indication 
for an extended population with a scale height of ~ 5 kpc. 

It is therefore natural to carry out a similar study for 
RR Lyrae stars since their nature is intermediate to RHB 
and sdB stars. This has become reliably possible using 
the recent magnitude limited {V < 12.5 mag) complete 
catalogue of RR Lyrae stars (Maintz 2005), a catalogue 
free of all kinds of misidentification. 

Sect.El describes the basic data and in Sect. 01 we 
present the results of our calculations. We analyse the ve- 
locities in Sect. 21 as well as correlations between velocity 
and other parameters like metallicity. The implications of 
our findings are discussed in Sect. Inland a comparison with 
results from other studies in Sect.El 

2. Data 

2.1. The sample 

The RRLyrs used for the present investigation are taken 
from the well documented sample established by Maintz 
(2005). She checked the identification of RRLyrs of the 
GSC with V < 12.5 mag against finding charts and with 
new observations. 

For the present investigation we restricted ourselves 
to those RRLyrs from her list which were observed in 
the Hipparcos mission (see The Hipparcos and Tycho 
Catalogues, ESA 1997; Hip and the Tycho 2 catalogue, 
H0g et al. 2000 Tyc2). Our sample thus has the best po- 
sitions and proper motions available. Furthermore, radial 
velocities have to be known. Of the 104 RRLyrs found 
in the Tyc2 we used 62, the other 42 having no reliable 
Vrad- From the Hip we could use 154 of the 182 RRLyrs 
for the same reason. Our data for radial velocity as well as 
distance from the galactic center came from sources in the 
literature (Beers et al. 2000, Fernley et al. 1998, Layden 
1994). We note that these authors used slightly different 
methods to come to stellar distances. We thus have inves- 
tigated 217 RRLyrs for their kinematic behaviour. Our 
sample has no kinematical bias. The limits are the limit- 
ing magnitude of the Hip and Tyc2 catalogues, combined 
with the availability of the Urad data. 

The present location of the RRLyrs of our sample in 
the X, Y and in the X, Z plane is shown in Fig.^ with 
hexagons. Evidently, RRLyrs have been searched for and 
found foremost at higher galactic latitudes, avoiding the 
crowded region of the disk. Furthermore, since our sample 
is magnitude limited, those stars included while suffer- 
ing from interstellar extinction have a smaller distance for 
their apparent magnitude. 

The RRLyrs being located now in the solar neigh- 
bourhood spread with time to nearly any galactic location 
(see Fig.n triangles). Many must have been at almost any 
location, both at very small and very large galactocen- 
tric radius and at small and large z-distance. This means, 
that our sample of 217 RRLyrs contains stars one can la- 
bel as part of a halo, of a thick disk, and of a thin disk 
population (not claiming these are sharp population dis- 
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Fig. 1. Distribution of the RRLyrs of our sample today 
(hexagons) and 100 Myr ago (triangles), equivalent to 
about half a galactic revolution. The full circle shows the 
position of the Sun 100 Myr ago. At left the distribution 
projected on the Galactic plane, at right perpendicular to 
the plane (along the X-axis). The circle in the left panel 
shows the 8.5 kpc radius. 



tinctions). The last category has kinematics keeping them 
confined to the galactic plane and thus they may be under- 
represented in our sample. Among the 217 are, however, 
several RRLyrs with solar abundances (see Sect. l2.2| and 
Sect. l3.2jl . suggesting that the lack of stars with thin disk 
kinematics need not be excessive. 



2.2. Metallicity of the RRLyrs 

The metallicity of our stars is taken from the litera- 
ture. We use data given by Beers et al. (2000), Fernley 
et al. (1998) and Layden (1994). As shown in Fig.El the 
metallicity of RRLyrs ranges from solar values to very 
low values of —2.6 dex. A broad peak in Fig.[51 shows 
that most of the stars of our sample have metallicities 
of —2 < [Fc/H] < —1. A small tail represents stars with 
still lower metallicity. These abundances indicate, follow- 
ing the usual arguments, that most of the stars are of pop- 
ulation II. But there are RRLyrs with higher metallicity 
as well: 46 stars have [Fe/H] > — 1. We have in our sam- 
ple 27 RRLyrs with a metallicity of [Fe/H] > —0.7 and 
11 stars even have [Fe/H] > —0.3. This suggests that the 
lack of relatively young RR Lyrs (likely having kinematics 
keeping them near the plane) cannot be excessive. 

3. Kinematics and orbits 

We have analysed the distribution of various kinematical 
aspects of our sample. To get the information for the over- 
all kinematics of the RR Lyrs and to attempt to come to 
a separation of disk and halo we calculated their orbits. 

3.1. Calculating velocities and orbits 

The observational parameters a, (5, d, /ia, ^5, Wrad were 
transformed into X, Y, Z, U, V, W (for details see de Boer 
et al. 1997b). Also the orbital velocity projected to the 
galactic plane, 6, and the velocity towards the galactic 
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Fig. 3. Meridional plots of some of the calculated RR Lyrae stars orbits. The orbits are calculated over a times pan 
of 10 Gyr. The top row shows orbits of thin disk star candidates. The middle row shows the orbit of RR Lyrae itself 
(on the left) and five more stars with elliptic orbits in/near the plane. The bottom row shows orbits of stars moving 
far out into the halo. For the orbit of the Sun see de Boer et al. (1997b, their figure 1). 
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Fig. 2. Histogram of the distribution of our stars in 
[Fe/H] . RR Lyrs cover a range from solar like values to very 
metal poor. Most stars show metallicity —2 < [Fe/H] < 
— 1. Note that several stars have metallicities near the so- 
lar value. 



centre, $, were calculated. Orbits were calculated over a 
time span of 10 Gyr^ taking steps of 1 Myr. Following de 
Boer et al. (1997b) the eccentricity, ecc, of the orbit and 



^ This long time allows to show more clearly the shape of the 
orbits. 



the apo- and perigalactic distances, i?a and _Rp, were cal 
culated as well as the maximum z-distance reached, z,„ax 
and its galactic radial distance normalised z-extent, nze, 

i?a — ^1 



ecc : 



and 



(1) 



For the values for our stars see Table ^ available in elec- 
tronic form. 

3.2. The orbits and the orbit parameters 

The meridional sections of the orbits of RRLyrs (see 
Fig.O show very different shapes. About half of the stars 
of our sample have perigalactic distances i?p < 3 kpc. 65 
RRLyrs reach a perigalactic distance i?p < Ikpc. Three 
stars (DRAnd, AMVir and AVVir) reach perigalactic 
distances less than 0.1 kpc. The other extreme is repre- 
sented by 5 RRLyrs (CI And, ALCMi, TWLyn, ARPer 
and SS Psc) which have perigalactic distances more than 
9 kpc and always stay beyond the solar circle. 

Of the 217 RRLyrs 63 have boxy orbits and stay close 
to the plane. A subset of 29 makes only small excursions 
in vj (Fig. 121 top row), 34 have planar but very eccentric 
orbits (Fig. O middle row). The orbit of RRLyr itself is 
not like that of a disk population star (even although its 
2max is only 0.21 kpc and its nze — 0.01), because its orbit 
has ecc = 0.8 with ~ 18.4kpc and Rp — 2.1 kpc. 

154 of the RRLyrs show orbits resembling those of the 
bottom row. The shape of these orbits is really chaotic, 
and shows movements perpendicular to the plane. These 
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Fig. 4. Histogram of the distribution of our stars in 8 
(00 = 220 kms~^). Some stars show circular velocities 
regarded to be characteristic for the thick and thin disk, 
respectively (see dashed lines). OlsRj ©thick disk and the 
velocity dispersions for thick and thin disk as found by 
Ojha et al. (1994) are indicated. Most RRLyrs are halo 
stars and the high peak represents the RRLyrs with ret- 
rograde motion. 



stars partly have orbits going to very small galactocentric 
distances and some of them reach very large apogalactic 
distances, with as an extreme i?p ~ 51kpc (for SXAqr). 

The statistics of the circular component Q of the ve- 
locity is shown in Fig.^J In our sample we have 87 RR Lyrs 
with retrograde orbits. We compared the characteristics of 
these stars with the prograde ones. Among the retrograde 
part of our sample there are no stars with a metallicity 
[Fe/H] > —0.9 (see Sect.0]upper left panel). Considering 
the retrograde group separately shows that the peak at 
high eccentricities in Fig.|Sl (left panel) is nearly com- 
pletely due to the retrograde RRLyrs, while the stars 
with prograde orbits show a flat distribution over the hole 
metallicity range. 

The eccentricities ecc of the orbits as well as the values 
of nze span a large range (see Fig.lSJ. 

The eccentricity of the orbit of the majority of our stars 
is ecc > 0.45. The distribution has an absolute maximum 
at ecc = 0.9 and a shallow local minimum at ecc ~ 0.6. 

The distribution of the normalised z-extents shows a 
maximum at low values nze < 0.2, a local minimum at 
nze — 0.4 followed by another peak at nze ~ 0.6. There 
are 109 RRLyrs with nze < 0.4 and only 45 RRLyrs 
reach nze > 1. A separation in stars with prograde and 
retrograde orbits does not change the appearance of the 
nze diagram. Both groups show high peaks at low values 
of nze and a long flat tail. 
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Fig. 5. Histograms showing the distributions of ecc (at 
left) and nze (at right) for the RR Lyrac stars of the 
sample. Note the flat distribution of ecc at low values and 
the peak at ecc ~ 0.9. Values of nze lie predominantly 
between and 1, with a local minimum at nze ~ 0.4 and 
a long tail up to nze — 3. 

3.3. Kinematic parameters and population type 

The shape of a stellar orbit gives an indication of the pop- 
ulation a star belongs to. Disk-population stars have orbits 
similar to those of the younger galactic stars, being rather 
circular and staying close to the plane. These are the ones 
most likely having been born in the disk. It is, from our 
data, not easy to discriminate between so-called thin- and 
thick-disk stars, a distinction which is not well defined 
anyway. The halo population consists of stars whose or- 
bits are not akin to disk-like galactic rotation. They most 
likely have been born outside the disk of the Milky Way. 

Using the above criteria (which were largely set by the 
distribution of the parameters of our stars as discussed 
above using Figs. 0] and EJ, one can now roughly sort the 
stars according to populations. We define halo stars as 
those having Q < 100 kms^^ or ecc > 0.4 or nze > 0.4. 
Each one of the criteria alone suffices to make the star a 
non-disk born star, and thus a star belonging to the halo 
population. 

Our statistics shows that out of the total of 217 stars, 
140 stars have Q < 100 km s^^, 153 stars have ecc > 0.4, 
and 108 stars have nze > 0.45, while 163 stars fulfil at 
least one of the criteria. Or, there are (kinematically) just 
54 classic disk stars in our sample. 

3.4. z-probability plot and scale height 

The sum of the orbits was used to derive the associated 
probability distribution of the z-distances of our RR Lyr 
star sample. This distribution is equivalent to the statisti- 
cal z-density gradient of the sample (for further details see 
de Boer et al. 1997b). One may fit the distribution with 
an exponential 



lniV(z)=lniVo--, 



(2) 



with iVo being the density in the local galactic plane and 
h the scale height. The slope of the log vs z distribution 
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Fig. 6. The probability of the z-distance distribution of 
aU stars of the sample, obtained by summing the individ- 
ual N{z) statistics for time steps of 1 Myr. The dashed 
lines show the two-component fit of the logarithmic dis- 
tribution based on the scale heights /i = 1.3 ± 0.1 kpc for 
the "disk component" and h = 4.5 ± 0.3 kpc for the "halo 
component" . 



gives the scale height. Fig.El shows our z-distance statis- 
tics. It turns out the left and right side are not identical 
(which is statistically not to be expected with the sam- 
ple size we have) but the differences are small. We fitted 
the left and right side separately. The left side gives scale 
heights of 1.25 ± 0.1 kpc for the steep central part and 
4.47 ± 0.3 kpc for the flanks. On the right side the val- 
ues found are 1.30 ± 0.1 kpc and 4.60 ± 0.3kpc, respec- 
tively. We arrived at these values after several tries, e.g., 
by changing the bin size and, because of small number 
statistics there, the limits of z to be included. For the fit 
of the central steep part we used a fitting interval of [0.3, 
2.9] and for the flanks [5, 33]. Given the similarity of the 
results we finally give one overall fit with 1.3 ± 0.1 kpc for 
the tight distribution and 4.5 ± 0.3 kpc for the extended 
distribution. The ratio of the mid-plane densities of the 
distributions is iV(0)haio/A^(0)disk = 0.16. 

The component with a steep z-distribution can be as- 
sociated with the disk population, both the thin and thick 
disk. We tried a 3 component fit to separate thin and thick 
disk but did not arrive at sensible results. It is, from this 
data, not possible to separate thin and thick disk. 

The analysis of the z-distribution statistics indicates 
the presence of (perhaps at least) two populations. The 
disk population consists of stars which stay near the plane, 
while the halo population ventures to larger z. In order to 
explore a possible sample bias in relation with the scale 
height of the halo population we have made two tests with 
subsamples, as follows. 

First, we selected those stars which at present are 
within the disk, i.e., those stars with Zpresent < 0.5 kpc. 
We made the z-statistics of this sample, as well as of the 
complement, and compared the two. The large scale z- 
distribution of the Zprcscnt < 0.5 kpc group shows a pro- 



nounced peak near z = kpc (which is of course what 
one expects). The z-distribution of the Zproscnt > 0.5 kpc 
group resembles the distribution of the whole sample but 
the mid-plane peak is missing. 

The flanks of the z-statistics toward large z of the two 
subsamples are very similar, indicating there is no sample 
selection in this sense, although the spatial distribution of 
the full sample is quite asymmetric (see Fig.P). 

In the subsample of the Zprescnt < 0.5 kpc group we find 
a slight separation of the thin and thick disk component 
leading to a scale height of 0.38 ± 0.03 kpc for the thin 
disk. This value is only slightly higher than that given by 
Chen et al. (2001). We note that our Zpj-cscnt < 0.5 kpc 
sample contains 63 stars. Still, the separation of thin and 
thick disk is not well defined in our sample. 

Second, we made a subdivision based on Zmax < 0.5 
kpc. Here the small z subsample is by definition restricted 
to stars staying close to the plane, while the Zmax > 0.5 kpc 
subsample must contain predominantly halo stars. The 
small z subsample (see Tab.|31l shows in a InA^ versus z 
diagram only a small high peak near z — kpc (by defini- 
tion). This is further discussed in Sect.El 

The large z subsample shows a z-statistics quite simi- 
lar to that of the full sample. It exhibits the scale height 
of 4.6 ± 0.3 kpc for the halo component, depending on the 
choice of bin size and the limits of z to be included. That 
value is hardly different from the halo scale height ob- 
tained from the full sample. Even the disk component is 
present with the same values of scale height as found for 
our full sample of RRLyrs, which is not surprising finding 
a scale height of the disk component of 1.3 kpc. The pres- 
ence of a disk component in this sample shows that the 
disk population is dominated by stars with orbits reaching 
beyond z = 0.5 kpc. 

We thus conclude that there is no significant sample 
bias in our RRLyr star sample. 

Altmann et al. (2004), dealing with sdB stars found 
with the same method a scale height of 0.9 kpc for the 
disk component and a halo component with a scale height 
of 7 kpc and Kaempf et al. (2005), considering RHB stars, 
give a scale height of 5 kpc for the halo component. 

4. Sample statistics of 6 with ecc, nze, [Fe/H] 

4.1. Circular velocity G 

The RR Lyrs of our sample have a mean circular velocity 
of = 47 km s^ ^ , but the stars cover a wide range (see also 
Tabled with +277 and —321 kms~^ as extreme velocities. 

Looking at the histogram of 8 (Fig.^J the most pro- 
nounced peak in the distribution is the one near kms^^, 
marking the group which in a broad sense does not par- 
ticipate in galactic rotation. 

The Toomre diagram (Fig.O allows to see the high 
asymmetry of our sample including many stars with pe- 
culiar velocities Wpoc > 100 km s^^. This asymmetry is the 
so called asymmetric drift (for more on that see, e.g., the 
review by Majewski 1993). RRLyrs with 9 < 6lsr domi- 
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Table 1. Mean values of V, W, 6, $, angular momentum, eccentricity and nze with dispersions for (subsets of) the 
sample. 



Subsample 


iV 


u 




V 




W 




9 




<1> 




J, 

J z 




ecc 


^ecc 


TLZE 
















[km s 


J 










[kpc km 


° J 










all 


217 




158 


+47 


125 


— o 


85 


+47 


123 


+5 


159 


+395 


1027 


0.61 


0.30 


0.63 


0.79 


rpe/Hl > —0 5 


20 


_2 


48 


+206 


31 


_3 


38 


+204 


32 


+3 


53 


+1718 


356 


0.19 


0.13 


0.09 


0.11 


Fe/H > -1 

L / J 


46 





82 


+165 


93 


-3 


53 


+168 


90 


+4 


82 


+1436 


810 


0.32 


0.28 


0.32 


0.67 


[Fe/H] < -1 


169 


-4 


173 


+15 


113 


+5 


92 


+ 14 


110 


+6 


175 


+113 


890 


0.69 


0.25 


0.71 


0.79 


[Fe/H] < -1.5 


86 


-25 


177 


-13 


114 


+5 


97 


-13 


108 


+27 


181 


-105 


880 


0.72 


0.23 


0.67 


0.61 


ecc < 0.45" 


58 


-15 


54 


+195 


40 


-19 


64 


+196 


40 


+19 


51 


+1644 


395 


0.22 


0.12 


0.19 


0.35 


ecc > 0.45 


149 


+ 1 


186 


-2 


92 


+11 


91 


-3 


88 





188 


-20 


722 


0.79 


0.17 


0.81 


0.85 


P < 0.35 d 


20 


-37 


133 


+44 


116 


+33 


100 


+44 


115 


+39 


134 


+387 


990 


0.59 


0.30 


0.64 


0.49 


0.35< P <0.55d 


103 


+2 


164 


+68 


137 


+2 


76 


+68 


135 


+1 


165 


+583 


1131 


0.57 


0.32 


0.53 


0.75 


P > 0.55 d 


91 


-1 


157 


+22 


108 


-1 


92 


+21 


106 


+3 


159 


+170 


861 


0.68 


0.26 


0.74 


0.87 



a) without seven stars having retrograde orbits 
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Fig. 7. Toomre diagram of Q versus Wporp, the ve- 
locity perpendicular to the Galactic plane, showing 
the asymmetry in kinematics of the RRLyrs. The 
dashed circles denote the peculiar velocity Vpcc = 
^/^^'^^W^~nQ~QhSR)^ in kms~^. The open circle 
shows the parameters of the Sun. 
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nate our sample and the dispersion of data points is large. 
Only few stars with Q > Qlsr are present. This is clearly 
correlated with the high eccentricity of the orbits of most 
of the RRLyrs as discussed in Sect. 13. 21 

The Bottlinger diagrams (see Fig.|SJ) show the asym- 
metry of data points as well. At low values of $ and W the 
RRLyrs do not concentrate near 0lsr but have a great 
variety of kinematics. The kinematic dispersion of W ex- 
ceeds that of $. The distribution of data points appears 
to be similar in the two panels. It seems to be symmetric 
in respect to both the $ = axis (upper panel) and the 
W = axis (lower panel). 



Fig. 8. Bottlinger and O ~ W diagrams show the orbital 
velocity of the RR Lyrae stars plotted against the velocity 
towards the Galactic centre ($, top) and perpendicular to 
the Galactic plane {W, bottom). A circle indicates the 
Solar values. 



4.2. Discussion of sample metallicities in relation with 
orbits and periods 

Given the orbit parameters, it is worthwhile to relate them 
with other parameters of RRLyrs. For that, we have di- 
vided our sample in subgroups in various ways (Table^. 
We used non kinematic selection criteria as metallicity and 
period but the kinematic criterion eccentricity as well. We 
choose as cuts [Fe/H] > —0.5, > —1, < —1, and < —1.5, 
for ecc < 0.45 and > 0.45, and for period P < 0.35 d, 
0.35 < P < 0.55 d and P > 0.55 d. 
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Table 2. [Fe/H] and Period with their dispersions for the 
217 RR Lyrae stars of the sample. 



Subsample 


TV 


[Fe/H] 


0"[Fe/Hl 


Period 
d 


fPoriod 


all 


217 


-1.33 


0.51 


0.51 


0.12 


[Fe/H] > -0.5 


20 


-0.28 


0.16 


0.44 


0.06 


[Fe/H] > -1 


46 


-0.55 


0.30 


0.44 


0.09 


[Fe/H] < -1 


169 


-1.54 


0.32 


0.53 


0.11 


[Fe/H] < -1.5 


86 


-1.80 


0.23 


0.56 


0.12 


ecc < 0.45" 


58 


-0.90 


0.57 


0.47 


0.09 


ecc > 0.45 


149 


-1.48 


0.36 


0.53 


0.11 


P < 0.35 d 


20 


-1.39 


0.48 


0.29 


0.07 


0.35< P <0.55d 


103 


-1.11 


0.51 


0.47 


0.05 


P > 0.55 d 


91 


-1.58 


0.40 


0.62 


0.06 



a) not including seven stars having retrograde orbits 



Metallicity is not a kinematic parameter. Yet, metal- 
licity gives no clear division between disk and halo stars. 
Stars with [Fe/H] > —1 may be members of the disk 
group in view of their Q (in particular those with [Fe/H] > 
—0.5). But the velocity dispersion ere — 90kms~^ is larger 
than expected for a pure disk population. So some of the 
RRLyrs with [Fe/H] > — 1 can be regarded to be halo 
stars. 

The more metal-poor RRLyrs with [Fe/H]< —1 are 
candidates for halo stars. This is corroborated by the re- 
sult of 8 = 14kms^^, (Te = llOkms^^, ecc ~ 0.69 and 
nze — 0.71 for this group. The velocity values of U, V 
and U for the metal poor part of our sample are similar 
to several published halo samples (see Martin & Morrison 
(1998), their Table 5) even if in their works the cuts for 
halo stars are set at more metal poor levels. Therefore we 
give the results for RRLyrs with [Fe/H]< —1.5, too. For 
this group our results shows to be in good agreement with 
results of previous investigations (see Martin & Morrison 
1998, Chen 1999). This metal-poor group of RRLyrs has 
Q = — 13kms~^, i.e., a mean retrograde rotation with 
respect to the LSR^. 

Orbit eccentricity supposedly divides between 
younger and older stars (more gravitational interactions 
of disk stars lead to larger deviations from the originally 
circular orbits). The subdivision at ecc = 0.45 shows 
the expected correlations with kinematic parameters (see 
Tabled notably Iz). This subdivision does not correlate 
in a pronounced way with metallicity. 

Note that we did not include 7 stars with retrograde 
orbits in the subgroup with ecc < 0.45 (Tablc^). These 
7 have a mean circular velocity Q = — 165kms~^ with 
dispersion of ce = 56kms~^. Also Chiba & Yoshii (1997) 
report the presence of halo stars with low eccentricity or- 
bits. Both results show that low eccentricity does not al- 
ways mean that a star is a disk member. The 58 RRLyrs 
are good candidates for disk stars (thin and thick disk). 



-200 200 

e I kins ' I 



0.2 0-4 0,6 0,8 1 
ecc 



1 



• * 



hi-fi:.: . 



^ Note that results for Q as given in the literature are some- 
times based on — 0, sometimes on Qps = 220 kms~^ 



10 20 30 0,2 0,4 0,6 0,8 

Z„. Per [d] 

Fig. 9. The plot shows metallicity [Fe/H] versus circu- 
lar velocity 8 (upper left panel), ecc (upper right), nze 
(lower left) and period (lower right panel). RRLyrs with 
[Fe/H] > — 1 show more disk-like orbits, but more metal- 
poor stars present both, disk-like orbits and halo-like or- 
bits having smaller 8, large eccentricity ecc and large nor- 
malised z-distance nze. With [Fe/H] versus period (lower 
right panel), note that there are 4 groups of stars. For the 
discussion see Sect. 14. 21 



Pulsation periods are related with the mass of the 
RRLyr, which is likely related with the development his- 
tory of the star. 

The groups with pulsation period P < 0.35 d and 
P > 0.55 d show both lower abundances than the large 
third group of periods 0.35 < P < 0.55 d. It has the largest 
abundance and the largest dispersion crpc/jj] . It is not a 
uniform subsample of RRLyrs. That this group is so dif- 
ferent is also visible in Tab.^ The dispersion in U, V, 
8, $, Iz and ecc is larger in this subsample than in the 
two other period groups. This fact can be noted as well in 
Fig. El (lower right panel). 

The group with periods of < 0.35 d contains the RRc 
stars with abundances of — 2 < [Fc/H] < — 1. The metallic- 
ity of the second period group ranges from solar to [Fe/H] 
of ~ —2. These are RRab stars. Note that the period in- 
creases in this group toward decreasing metallicity. Most 
of the RRLyrs of our sample are in this group. It com- 
bines RRLyrs with a large dispersion in metallicity and 
various orbit parameters, as can be recognised in Tab. El 
probably representing different ages. The third group hav- 
ing periods ~ 0.6 d are RRab stars as well but with lower 
abundances. This group is the field RRLyr equivalent to 
the RRLyrs in Oosterhoff II clusters. This last group is a 
very small subsample of RR Lyrs with relative long peri- 
ods and metallicities [Fe/H] < —1.6. 
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All parameters are considered in Table[21 RRLyrs 
with ecc > 0.45 have lower abundances ([Fe/H] = —1.48) 
than those with ecc < 0.45 ([Fe/H] = —0.90), as we ex- 
pected. Like in Tab.^we did not include 7 RRLyrs with 
retrograde orbits having ecc < 0.45. We did so because 
these stars do not fit into the group of low eccentricity 
stars which are candidate disk stars. This group of retro- 
grade RRLyrs has a low metallicity of [Fe/H] = —1.85. 

Fig.ini shows the data of Tables and 13 in graphi- 
cal form: metallicity [Fe/H] versus orbital velocity Q, ecc, 
nze, and period. Dividing the RRLyrs in metal-poor and 
metal-rich groups as Layden (1995), we also find that ve- 
locities for RRLyrs with [Fe/H] > — 1 reveal a distribution 
different than those with lower metallicities. High metal- 
licity stars with more disk-like circular velocities O clus- 
ter both at low ecc and nze. But the separation is not 
a strict one. Even some data points with very low abun- 
dance [Fe/H] < -1.5 show 9 ~ 200 kms~^, ecc < 0.4 and 
nze < 0.2. Altmann & de Boer (2000) found the same fact 
using a much smaller sample of RRLyrs. 

These findings are similar to those by Chiba & Yoshii 
(1998) who report that even stars with very low metal- 
licity may have "nearly circular orbits with ecc < 0.4". 
High metallicity does indicate disk-like kinematics, but 
low metallicity does all by itself not give sufficient infor- 
mation about the orbit shape. An eccentricity of < 0.4 
is alone not a sufficient criterion for population. Chiba & 
Beers (2000) find in their investigation of metal-poor stars 
"a remarkable discontinuity of the rotational properties of 
the Galaxy" at [Fe/H] ~ —1.7 (see also their figure 3). 
Our sample of RRLyrs does not confirm this. We find a 
continuous distribution in O, ecc and nze in the metal- 
licity range — 2 < [Fe/H] < —1. But we note that in a 
sample of 163 metal-poor halo stars chances are non-zero 
to find one with an orbit rather co-planar with the Milky 
Way disk. 

5. Discussion of results 

We calculated the orbits of the RRLyrs of our sample 
and derived the z-distribution statistics. This led to scale 
heights for a disk group of RRLyrs (/idisk — 1.3 kpc) and 
for a halo group (/ihaio — 4.5 kpc). However, analysing 
for each star the parameters Q, ecc, nze, and [Fe/H] it 
became clear that all parameters are needed to classify a 
given star as a genuine disk or halo star. Of our 217 stars, 
63 are truly disk like, 154 are halo like. 

In our sample we have 87 RRLyrs with retrograde or- 
bits. We compared the characteristics of these stars with 
the prograde ones. Among the retrograde part of our sam- 
ple there are no stars with metallicity [Fe/H] > —0.9 (see 
Fig.|n|upper left panel). These retrograde rotating RRLyrs 
are not separated from the rest of the sample in respect 
to z,nax or period, but their group values of ecc — 0.79, 
Cecc — 0.20, nze — 0.79 and anze = 0.79 are higher than 
the values for prograde halo candidates of about the same 
metallicity. We again conclude that most of the RRLyrs 
of our sample are halo stars. 



Table 3. 16 RRLyrs staying close to the plane 
(zmax<0.5kpc) 

Name [Fe/H] O Ra Rp Zmax ecc 

[kms^''"] [kpc] [kpc] [kpc] 
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Fig. 10. Plots ecc against [Fe/H] and Zmax- On the left 
side only RRLyrs with prograde orbits (8 > 0) are plot- 
ted, on the right side only those with retrograde orbits 
(6 < 0). The eccentricities of the prograde RRLyrs spread 
over the hole range of ecc while the retrograde ones show 
foremost higher eccentricities. 



A different way of looking at the data is given in 
Fig. El There [Fe/H] and Z^ax are plotted against ecc for 
the prograde and retrograde RRLyrs. Clearly, the contin- 
uous distribution indicated above is present for the pro- 
grade set, while the retrograde stars rather are metal poor 
and have orbits with predominantly larger eccentricity. In 
the retrograde set there are only stars with [Fc/H] < —0.8 
and ecc > 0.45. We conclude from this that the prograde 
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set does contain a good fraction of higher metaUicity low- 
eccentricity orbit stars, in short younger RRLyrs, while 
the retrograde set has only those with characteristics not 
like that of a younger galactic disk. 

We took a special look at the stars not reaching beyond 
^max = 0.5 kpc These stars are listed in Tab.|21 Among 
these stars are two with a much smaller metallicity than 
the rest and they turn out to have very eccentric orbits. 
These are RR Lyr itself, and V338 Pup whose orbit is ret- 
rograde. Thus 2 of 16 stars staying close to the plane turn 
out not to have disk-like orbits and thus belong to the 
halo population. 

Chiba & Beers (2000) published a list of likely metal- 
weak thick disk stars including 3 RRLyrs using as cri- 
terium [Fe/H] < -1, ecc < 0.2 and 183 < < 225 kms"!. 
In our sample we found 7 such RRLyrs: the three found 
by Chiba & Beers (2000) plus TZAqr, RXEri, V690Sco 
and AT Ser. We confirm that there are metal-weak stars 
with classical disk orbits, which may be part of a metal- 
weak thick disk or in two cases (V494 Sco and V690 Sco) 
even thin disk star candidates having a z,iiax<0.5kpc. 

6. Halo distribution of RR Lyrae stars 

The spatial distribution of RRLyrs and other HB-like 
stars has been investigated several times based on di- 
verse samples. In a pioneering study. Plant (1965) de- 
rived a scale height of RR Lyr stars in the solar vicin- 
ity of ~2.3 kpc. In his review Majewski (1993) reports 
scale heights for the "intermediate population 11" , being 
~ 1.5 kpc for RRLyrs. In almost all cases, the studies 
used photometric distances perhaps together with radial 
velocity information. Proper motions have thus far been 
included in limited cases only. 

The vertical scale height derived in this paper from the 
orbit statistics can be used for a comparison with results 
from other studies. 

Kinman et al. (1994) used BHB stars (with colours 
between the blue edge of the RRLyr strip and the inter- 
section of the HB with the main sequence, so rather HBA 
stars) to derive the structure of the halo from statistically 
complete samples. They assume a (spherical) radial dis- 
tribution of the form 

P - P0i?gal ' (3) 

and find that the distribution of their HBA stars agrees 
with that. We have fitted our z-distribution to the baro- 
metric formula (Eq.[21) and arrived at a scale height of 
~ 5 kpc (see Fig.|Hl). Transforming the run of densities of 
a radial distribution into an exponential vertical distribu- 
tion at the location of the Sun {-w ~ 8.5 kpc) we find that 
the radial and the barometric relation (with h — b kpc) 
give a run of densities N{z) which, for 5 < z < 20 kpc, de- 
viate from each other by only 20% (the barometric formula 
being a bit steeper). Thus our result from RRLyr star or- 
bits shows a spatial distribution similar to the distribution 
of HBA stars from samples by Kinman et al. (1994). 



The mean space density of RR Lyrs in the solar neigh- 
bourhood has been derived by Amrose & Mckay (2001). 
They arrive at different values for different solar distances 
but their overall result is fairly summarised by p = 5 ± 1 
stars kpc~^. They also derive a vertical distribution and 
give fits with the barometric formula allowing scale heights 
between 0.75 and 1.8 kpc (using a single exponential fit). 
However, from their Figure 5 it is clear that the RRLyrs 
show a steeper distribution at low z with /i ~ 1 kpc, simi- 
lar to our result, and a flatter one further out of /i > 2 kpc. 

Our data cannot be used to find an absolute mid- 
plane density. However, one could use the ratio of 
A^(0)haio/A^(0)diskRRLyrs = 0.16 together with the Amrose 
& Mckay value of p = 5 ± 1 stars kpc""^. Furthermore, the 
midplane ratio may be combined with the scale heights de- 
rived for the disk and halo components to get the ratio of 
all RRLyrs in these components. The procedure is similar 
to that used by de Boer (2004) for sdB stars. We so find 
that the Milky Way has about half as many halo RR Lyrs 
than disk RRLyrs. Since the midplane density ratio prob- 
ably is still affected by systematic effects and the scale 
heights are, among various studies, only just consistent, 
the indicated total RRLyr star ratio has to be regarded 
as a crude value. 

The data from the Sloan Digital Sky Survey (SDSS) 
revealed numerous HB stars also at very large distances 
from the Milky Way disk (see, e.g., Yanny et al. 2000, 
Sirko et al. 2004). Vivas et al. (2001), using data from the 
QUEST survey for RR Lyrs, find an over-density of RR Lyrs 
at 50 kpc from the Sun. They surmise this over-density is 
part of the tidal stream associated with the Sagittarius 
dwarf spheroidal galaxy. From these new survey results 
we conclude that the halo is, on a large scale, apparently 
rather inhomogeneous so that RR Lyr distribution values 
will show to vary in different directions and at different 
distances. 

7. Conclusions 

RRLyrs exist in the Milky Way as part of the disk and 
part of the halo. We separated the halo component from 
disk candidates using the shape of their orbits. Of 217 
RRLyrs, 154 or 71 % are member of the halo population. 
Therein we find a retrograde rotating halo component 
of 87 RRLyrs or more than half of our halo group. 
The (vertical) spatial distribution found with our orbit 
statistics method is, out to z =20 kpc, similar to that 
found in other studies using classical methods. We show 
that there are RRLyrs staying in the disk but which are 
member of the halo population. 
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Table 4. Electronic Table; Positions, velocities (given 
in the Galactic euclidic system XYZUVW and 9,$) 
and orbital data (iia, -Rp, -^max, normalised z-extent, nze 
and eccentricy ecc) of the stars of our sample. This ta- 
ble is only available in electronic form at the CDS via 
anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat7J/A-l-A/ 



Name XYZUVW^OR^^Rp z^ax nze ecc 

[kpc] [kms""""] [kpc] — — 

Online-Data are available under 
http; / /www.as tro.uni-bonn.de/~gmaintz| 



